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Abstract 

The first goal of Vibration- Transit (V-T) theory was to construct a tractable approximate Hamil- 
tonian from which the equilibrium thermodynamic properties of monatomic liquids can be calcu- 
lated. The Hamiltonian for vibrations in an infinitely extended harmonic random valley, together 
with the universal multiplicity of such valleys, gives an accurate first-principles account of the 
measured thermodynamic properties of the elemental liquids at melt. In the present paper, V-T 
theory is extended to non-equilibrium properties, through an application to the dynamic structure 
factor S{q,uj). It was previously shown that the vibrational contribution alone accurately accounts 
for the Brillouin peak dispersion curve for liquid sodium, as compared both with MD calculations 
and inelastic x-ray scattering data. Here it is argued that the major effects of transits will be 
to disrupt correlations within the normal mode vibrational motion, and to provide an additional 
source of inelastic scattering. We construct a parameterized model for these effects, and show 
that it is capable of fitting MD results for S{q,uj) in liquid sodium. A small discrepancy between 
model and MD at large q is attributed to multimode vibrational scattering. In comparison, mode 
coupling theory formulates S{q,uj) in terms of processes through which density fluctuations decay. 
While mode coupling theory is also capable of modeling S{q,uj) very well, V-T theory is the more 
universal since it expresses all statistical averages, thermodynamic functions and time correlation 
functions alike, in terms of the same motional constituents, vibrations and transits. 

PACS numbers: 05.20.Jj, 63.50.+X, 61.20.Lc, 61.12.Bt 
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I. INTRODUCTION 



Vibration- Transit (V-T) theory is a Hamiltonian formulation of the dynamics of 
monatomic hquids. It is based on the idea that a hquid system moves on a potential 
energy surface making jumps between valleys, that these jumps are approximately instan- 
taneous, and that the dominant majority of visited valleys are all random in structure and 
are equivalent in energy and vibrational properties. The zeroth order approximation to the 
Hamiltonian expresses the liquid motion in terms of normal mode vibrations in a single 
infinitely-extended harmonic random valley, and can be explicitly calculated from first prin- 
ciples for actual systems. It is now well known that this vibrational motion gives a very 

Biod account of the equilibrium thermodynamics of monatomic elemental liquids at melt 
This result is conceptually fundamental, not only because it supports the po- 
tential energy landscape picture of liquid dynamics, but also because it validates the basic 
assumptions of V-T theory. Moreover, it was obtained by V-T theory without adjustable 
parameters, a result that no other tractable theory has achieved. Following this success, in 
the present work we look into a deeper level of the dynamical behavior in liquids, namely 
its nonequilibrium properties, and apply V-T theory to time correlation functions, which 
express nonequilibrium properties in linear response theory Again without adjustable 
parameters, the vibrational contribution to any time correlation function can be calculated 
from the zeroth order Hamiltonian. Once more the vibrational contribution turns out to 
play a central role and precisely gives the location of the Brillouin peak in the inelastic 
scattering data for liquid sodium 5|. However, the width of the Brillouin peak is larger than 
the vibrational width alone. This broadening of the Brillouin peak results from the transit 
contribution, for which an explicit evaluation is not yet available. The purpose of this paper 
is to construct and test a model for the contribution of transits to inelastic scattering. The 
model is shown to be very successful for our case study, and thus provides a new description 
of the scattering process. 

The dynamics of liquids and supercooled liquids, studied with the aid of MD calculations 
for glass forming systems, is a currently active research field. We have presented extensive 
comparison of V-T theory with a broad range of potential-energy-landscape theories . In 
a paper of particular relevance here, Mazzacurati, Ruocco, and Sampoli (see also Qj) 
have shown that a vibrational analysis is in excellent agreement with MD calculations of 
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S{q,uj) for a Lennard- Jones glass. Beyond this result, extensive theoretical analysis of the 
complete atomic motion is required before a vibrational contribution can be incorporated 
into a theory of liquid dynamics. This analysis contitues the foundation of V-T theory 
and provides the following stipulations: (a) potential energy valleys used in liquid 
theory must be random valleys, and not some other symmetry; (b) because all random valleys 
of a given system are equivalent in vibrational properties, the liquid vibrational contribution 
can be calculated from a single random valley; (c) the representative random valley has 
to be extended to infinity so that the vibrational statistical averages are defined; (d) for 
thermodynamic functions, corrections for anharmonicity and valley- valley intersections must 
be recognized; and (e) for time correlation functions, the vibrational motion has to be 
supplemented with transits in the liquid. These stipulations are crucial to the present 
theoretical development. 

Starting from the exact first-principles vibrational contribution, our model for transits 
in dynamic response theory is constructed in Sec. II. In Sec. Ill, the model parameters are 
adjusted to achieve agreement with MD calculations for a system representing liquid sodium 
at melt. The quality of the fitting is discussed, as well as the interpretation of the fitted 
parameters. For this example of dynamic response in monatomic liquids, we present in 
Sec. IV a detailed comparison of our V-T theory with mode coupling theory, which is at 
present the most successful in accounting for time correlation functions in liquids. Our 
conclusions are summarized in Sec. V, and the unifying nature of V-T theory is noted. 

II. CONSTRUCTION OF THE TRANSIT MODEL 

Let us consider a system of atoms in a cubic box with the motion governed by periodic 
boundary conditions. The position of atom K at time t is rK{t), K = 1, . . . ,N. The density 
autocorrelation function is 



where (...) represents a thermal average over the motion, plus an average over the star 
of q, which converts the right side to a function of q for finite systems. In V-T theory of 
the liquid state, the motion consists of normal mode vibrations within a single extended 
(harmonic) random valley, plus transits between valleys. We shall neglect anharmonicity. 




(1) 
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and will consider classical motion so that position coordinates may be commuted at will. 

Let us neglect transits for the moment, and consider motion in a single random valley. It 
is convenient to write 

rK{t) = RK + UK{t), (2) 

where Rk is the equilibrium position and uxit) the displacement. The contribution to 
t) is the vibrational contribution, given by 

FUq, ^) = ^ E (e-^^-^^-W-^^^"))) , (3) 

KL 

where Rrl = Rk^Rl- The motional average is now a harmonic vibrational average (. . . )^, 
and Eq. (jS)) simplifies to (see e.g. 

FvUq,t) = ^( E^"^"'"""' e-^-('')e-'^-('^) [1 + {q ■ q ■ n^(0)), + ...]), (4) 

\ KL I q* 

where Wxiq) is the Debye- Waller factor for atom 

WK{q) = \{{q-UKf)^, (5) 

and where (. . . )g* is the q-star average. The series in brackets in Eq. is the expansion of 
an exponential. Since {q ■ Uxit) q ■ ul{0))i^ vanishes as t — > oo, the constant term in Eq. (jlj) 
is Fyib{q, oo), given by 



oo) = / 5^ cos(q ■ Rkl) e-^^^-^^e" 

\ KL 



(6) 



where cos(q ■ Rkl) appears because of the star average. To leading order in the expansion 
in Eq. (jH), the time dependence of F^ih{q,t) is contained in the function 

- oo) = V 5^cos(q ■ Rkl) e'^^^'^^e-'^^^'') {q ■ ukU) q ■ u^m^ \ . 

\ KL I q* 

(7) 

To evaluate the dynamic structure factor S{q,uj), the displacements are written as a 
sum over normal modes A, which have frequencies uox and eigenvector components Wkx-, for 
A = 1, . . . , 3A^. The result is joj 

S^ib{q,uj) = F^ib{q,oo)6{uj) + Sl]l{q,uj), (8) 
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where 



A 



/a(?) = -( >:e-'^-"-e-™q-^;,A ) • (10) 

The first term in Eq. (jHJ describes elastic scattering, while Slj]l{q,uj) describes inelastic 
scattering from the vibrational normal modes in the one-mode approximation. Multimode 
scattering will arise from the higher order terms in Eq. (j3)), and are neglected in the present 
work. It is understood that the three modes of uniform translation, for which ux = 0, are 
omitted from all statistical mechanics equations. 

We now allow for transits. When atom K is involved in a transit, both Rk and Uk change 
in a very short time, in such a way that Rk + remains continuous and differentiable 
in time. A detailed model of transits in the atomic trajectory may be found in Chisolm et 
al. 0,0. Here we seek a simpler approximation. If the time segments between transits 
involving atom K are denoted jk = 1,2, ... , then the position of atom K at time t is 
RKilxit)) + UKijKif),^)- F{q,t) for the liquid is then written, from Eq. ((H), 

Our numerical studies provide evidence, described below in connection with Eq. ()14|) . that 
transits can be approximately neglected in the displacements UK{'yK{t),t). We therefore 
make this approximation, and separately average the displacement terms in Eq. flllj) over 
harmonic vibrations. The result is Eq. with R^ replaced by Rxi^Kif))) 

Fu,{q,t) ^ l / ^e-''-[^-(^-W)-^-Wl e'^-^'^^e-'^-^'^) [1 + {q ■ UKit) q ■ ^^z.(0)), + . . . 

\ KL 

;i2) 

Our next step is to modify this equation so as to model the presence of transits in Rk- 

There are two ways in which transits contribute to Fiig{q,t). First, transits introduce a 
fluctuating phase in the complex exponential in Eq. p2|l . and this causes additional time 
decay through decorrelation along each atomic trajectory. We model this with a relaxation 
function of the form e~"*. Second, transits give rise to inelastic scattering, in addition to 
the vibrational mode scattering already present in Eq. ()12|) . and this increases the total 
scattering cross section. We model this with a multiplicative factor. 



The leading term in Eq. (fT^ gives rise to the hquid Rayleigh peak, and so is denoted 
FniQyt)- Without transits, FR{q,t) reduces to Fmb{q,oo), Eq. (0), so we model FR{q,t) as 

Fn{q, t) = C{q) F,,t,iq, oo) e-"^^'?) *. (13) 

This function decays to zero with increasing time, in accord with the liquid property 
Fiiq{q,t) — s> as t — s> oo. The relaxation rate ai{q) is expected to be around the mean 
single-atom transit rate. C{q) is positive, and greater than 1 because of the inelastic scat- 
tering associated with transits (notice the total scattering cross section is not affected by 
the factor e~"^*). 

The displacement-displacement correlation function in Eq. (fT^ gives rise to the Brillouin 
peak, and so is denoted FB^q, t). Without transits Fsiq, t) reduces to Fmb{q, t) — Fmb{q, oo), 
Eq. ((7j). Empirically, we have found that this vibrational contribution alone gives an ex- 
cellent account of the location of the Brillouin peak, and the total cross section within it, 
as compared with MD calculations and with experimental data for liquid sodium j^. This 
suggests keeping the vibrational contribution intact, as we did in going to Eq. (fT2)) . and also 
suggests that we model ^^(g, t) by 

FB{q,t) = [F„,(g,t) -F„,(g,oo)] e'^'^^^ K (14) 

Note Fmb{q, t) — Fmbig, oo) decays to zero with time, because of the decay of the vibrational 
correlation function in Eq. (jZj), and this decay gives the Brillouin peak its natural width j^. 
The right side of Eq. ()14p decays faster with time, hence broadens the Brillouin peak from 
its natural width, but leaves its total cross section unchanged. 

From the above equations, our model for the dynamic structure factor is 

Shg{q, uj) = Sn{q, uj) + Ssiq, t^), (15) 



c. / N C{q) ai F^ibjq, oo) 

SR{q,uj) = — Y\ ' (Id) 



7r(a;2 -|- af) 

A 



^2 «2 



The model has three adjustable parameters for each q, namely C{q),ai{q), and a2{q)- 



(17) 
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FIG. 1: Vibrational mode frequency distribution for a single random valley (1497 modes). 



We note the presence of a short-time error in Fiig{q, t). The correct short-time behavior is 
Fiiq{q, 0)+b . . , with known coefficient & ( J], Eq. (7.4.41)). The vibrational contribution, 
Eq. ((7j), has the correct limiting behavior, but the model functions in Eq. (fT^ and (fT^ are 
linear in time, since e"""* = 1 — at + ... . The linear term is important up to a time tc, which 
is very small, and beyond tc the time dependence of Fmb{(l,t) dominates. In our system 
we estimate the linear time dependence contributes to Suq^q.uj) only at frequencies above 
50ps~^, which is above the largest uj\ present. For reference, the normal mode frequency 
distribution g{uj) is shown in Fig. 1. 

To complete this Section, let us estimate the average rate v at which an individual 
atom is involved in a 

our general estimate for monatomic liquids at melt is z/ ~ {uo) /27r, where (uj) is the rms 
vibrational mode frequency. For liquid sodium at melt, this gives v ~ 2.5 ps^^. 



III. RESULTS AND DISCUSSION 



The system we study has = 500 atoms with an interatomic potential representing 
metallic sodium at the density of the liquid at melt. The potential gives an accurate account 
of the vibrational and thermodynamic properties of crystal and liquid phases, and a good 
account of self diffusion in the liquid (for summaries see Here we use the sodium 
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TABLE I: Quantities associated with our liquid S{q,uj) model. 



Q 


Q (ao^) 


Fyib{q, oo) 


F^l>{q,0)/F,,,{q,0) 


C{q) 


(0,0,1) 


0.12129 


0.0043 


1.00 


2.0 


(1,1,1) 


0.21009 


0.0033 


0.98 


2.3 


(1,1,2) 


0.29711 


0.0026 


0.94 


2.9 


(1,1,3) 


0.40229 


0.0021 


0.90 


4.0 


(0,1,4), (2,2,3) 


0.50011 


0.0021 


0.84 


5.1 


(0,3,5), (3,3,4) 


0.70726 


0.0044 


0.72 


8.3 


(0,3,6), (2,4,5) 


0.81367 


0.0110 


0.70 


9.4 



potential to see how well our transit model can be made to fit MD results for S{q,uj). We 
study g-values in the range from 0.12 ag^, the smallest allowed q for our system, up to 
0.81 aQ^, beyond which the Brillouin peak is poorly discernable. In comparison, the first 
peak in S{q) is at qm = 1-05 a^^. The model is evaluated from Eqs. P5II17|) for a single 
random valley, and the results show scatter due to the small system size. To reduce this 
scatter we used a graphically smoothed curve of F„jf,(g, oo) in Eq. (fTBj) for SR{q,uj), where 
the smoothed data are listed in Table I. In comparison, the MD results show little finite- iV 
scatter, since the MD system visits a very large number of random valleys during the decay 
time of F{q, t). 

The individual functions SR^q, u), Ssiq, ^^), their sum Suq^q, u), and SMoiq, ^^) are shown 
in Fig. 2 for a representative q. The slight decrease in SB{q,i^) at small q appears because 
our system has no vibrational modes with frequencies below 1.7 ps^^ (see Fig. 1). In the 
fitting process, we adjusted C and ai to get fits of the intercept at = 0, and of the slope 
in the steeply decreasing range at u ^ 1-3 ps~^, and we adjusted a2 to get an overall fit to 
the Brillouin peak. 

The fitted Sug^q, u) and Smd{(1, ^) are shown for the remaining g- values in Figs. 3-8. For 
each q, the shapes of the components SR{q,uj) and SB{q,uj) are qualitatively the same as 
those shown in Fig. 2. The overall fits of our model to MD data are very good for q from 
0.12 to 0.50 aQ ^, from Figs. 2-6. Except for the small undershoot of theory at the bottom 
of the dip between Rayleigh and Brillouin peaks, and the overshoot of theory in the high 
frequency tail, the model discrepancies can be attributed to scatter due to evaluation for 
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co(ps ) 



FIG. 2: S{q,uj) for q = 0.12129 ag^, from the model (solid line) and from MD (circles). The 
Rayleigh (broken line) and the Brillouin (dotted line) contributions to the model are also shown 
separately. 

only one random valley. But at q around 0.71 and 0.81 a^^, Figs. 7 and 8, the model cannot 
be made to fit the Brillouin peak quite so well as in the preceeding five figures. The figures 
suggest that our overall physical description is still correct at the two highest q, but a small 
correction needs to be addressed there. 
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FIG. 3: Same as Fig. 2 but for q = 0.21009 a^ 
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FIG. 4: Same as Fig. 2 but for q = 0.29711 ag ^ 
12 
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FIG. 5: Same as Fig. 2 but for q = 0.40029 aQ \ 

Experience with inelastic neutron scattering in crystals suggests that multimode scat- 
tering should become significant at temperatures close to melting and q beyond the first 
Brillouin zone boundary. This includes our system at g > 0.70 a^^ . We note that oo), 
Eq. is not affected by the one- mode approximation. Since 0) is the maximum 

magnitude of t), an estimate of the accuracy of the one- mode approximation is pro- 

vided by the ratio F^l^{q, 0)/Fyih{q, 0), where the numerator is the one-mode approximation 
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FIG. 7: Same as Fig. 2 but for q = 0.70726 a^^ . 



and the denominator is exact. The ratio is hsted for each q in Table I. We interpret the 
results as follows. Multimode scattering, not included in Eq. (fTTj) for Sb{<1,^^), is present 
in the MD data, is small for q < 0.50 a^^, but is responsible for the discrepancy between 
model and MD around the Brillouin peak in Figs. 7 and 8. 

Let us now consider the magnitude of the fitted parameters. The Rayleigh peak strength 
C{q), listed in Table I, is greater than 1 and increases steadily as q increases. This implies a 
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considerable cross section for inelastic transit scattering. The Rayleigh peak relaxation rate 
ai{q) is graphed in Fig. 9. From Eq. (fTB|) . ai{q) is the Rayleigh peak half width at half max. 
ai{q) extrapolates toward zero as g ^ 0, while ai{q) is roughly constant at g > 0.5 a,Q^. 
The Brillouin peak relaxation rate a2{q) is also graphed in Fig. 9, and appears to go to 
zero at q around 0.1 a.^^. Notice 02(5) does not measure the width of the Brillouin peak, 
but measures its width beyond the natural width. Hence the 0:2(1?) graph suggests that the 
liquid Brillouin peak has its natural (vibrational only) width at small 5', at g < 0.1 a,Q^ in 
the present work. Except where the relaxation rates approach zero at small q, a\ and 02 
are in the range 1-5 ps~^, in qualitative agreement with the mean transit rate v ~ 2.5 ps^^. 



IV. COMPARISON WITH MODE COUPLING THEORY 



The theory most successful to date in accounting for time correlation functions is mode 



coupling theory |11|. Detailed summaries of mode coupling theories of liquid dynamics are 
given by Boon and Yip [12], Hansen and McDonald ^, and Balucani and Zoppi Q. Mode 



coupling theory has been applied to the glass transition 



15, 



16|, and has been shown 



capable of .ational^tag the den.ty co.elat.on function, at te.pe.atu.e. . the vicinity of 



the glass transition 



17L llSL Il9j . The application for which we shall compare mode coupling 



and V-T theories is dynamic response of monatomic liquids at temperatures near and above 
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FIG. 9: Relaxation rates from our liquid model: a\ (filled circles) for the Rayleigh peak 
and ai (empty circles) for the Brillouin peak. 



melting. 

Mode coupling theory works with the generalized Langevin equation for F(g, t), and 
expresses the memory function in terms of processes through which density fluctuations 
decay [J, In the viscoelastic approximation, the memory function decays with 

a g-dependent relaxation time This approximation provides a good fit to the 

combined experimental data [2^ andMD data j^, 22| for the Brillouin peak dispersion curve 
in liquid Rb j23| (see also Fig. 9.2 of [4]). Going beyond the viscoelastic approximation, Bosse 
et al 0, 2^ constructed a self-consistent theory for the longitudinal and transverse current 
fluctuation spectra, each expressed in terms of relaxation kernels approximated by decay 
integrals which couple the longitudinal and transverse excitations. This theory is in good 
overall agreement with extensive neutron scattering data and MD calculations for Ar near 
its triple point 25j. The theory was developed further by Sjogren |2a, |22[, who separated 
the memory function into a binary collision part, approximated with a Gaussian ansatz, 
and a more collective tail represented by a mode coupling term. For liquid Rb, this theory 
ives an "almost quantitative" agreement with results from neutron scattering experiments 



2011 and MD calculations 



21 



221. More recently, inelastic x-ray scattering measurements 
Kave been done fo. the a.Ua„ ,neU>s L, Q and Na Q fi. These data have bee. 

analyzed by mode coupling theory, and the resulting fits to ^(g, tu) are excellent, both for 
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the experimental data and for MD calculations l29l |31 

A detailed comparison of the present study with the analysis of Scopigno et al. j29[ is of 
interest. They analyzed experimental data for liquid sodium at 390 K, for q in the range 
0.08-0.77 ao ^ Their memory function has three relaxation terms, one for coupling between 
thermal and density degrees of freedom, with no adjustable parameters, and two for true 
viscous processes, each having adjustable weight and relaxation time. One parameter is 
fixed by the total weight, so Scopigno et al. have effectively three g-dependent parameters 
to fit the shape of S{q,uj). In the present work, we calibrate V-T theory by comparison 
with MD data for liquid sodium at 395 K, for q in the range 0.12-0.81 ag ^. Our Rayleigh 
peak contribution has one weight parameter and one relaxation rate, and our Brillouin peak 
contribution has one relaxation rate, making three adjustable g-dependent parameters. Our 
Brillouin peak weight parameter was fixed at one because we had already confirmed that 
the vibrational contribution alone has the correct weight. We find that our ai is nearly the 
same as Scopigno et al.'s inverse relaxation time r^^ over the entire q range: the relative 
difference in magnitude averages 25%. Though we understand the reason, that in each case 
the parameter is determined by the width of the Rayleigh peak, the level of agreement is 
remarkable nevertheless. On the other hand, our 0^2 is much smaller than their r~^, as 
the ratio a2/T^^ varies from 0.01 at small q to 0.07 at large q. Again the reason is clear: 
while is determined by the Brillouin peak width, ^2 is determined by only the width 
beyond the natural width. The result illustrates the important point of comparison between 
V-T and mode coupling theories: the two methods are based on different decompositions 
of the physical processes involved. While mode coupling theory analyzes F{q,t) in terms of 
processes by which density fluctuations decay, V-T theory analyzes F{q,t) in terms of the 
two contributions to the total liquid motion, vibrations and transits. 



V. CONCLUSIONS 



By resolving the complete atomic motion into its constituents, vibrations and transits, V- 
T theory offers a unified theoretical formulation of equilibrium statistical mechanics averages, 
both of thermodynamic variables and of time correlation functions. The vibrations alone 
provide a tractable, accurate, parameter free formulation of thermodynamic properties of 
monatomic liquids [l|. The role of transits is merely to allow the liquid to visit the vast array 
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of random valleys, and hence achieve the full liquid entropy 

HQS 

. The same vibrations 

provide a tractable, parameter free contribution to time correlation functions, while the 
same transits, which are part of the equilibrium fluctuations, complete the formulation of 
time correlation functions. 

Our study of the dynamic structure factor exemplifies this unification. The vibrational 
contribution provides the nontrivial results in Eqs. (pirUj) for Smbiq, uj). Pure elastic scatter- 
ing is given by Fmb{q, oo) 6{lj), where Fmb{q, oo) is the positive long-time limit of Fmb{q,t). 
Inelastic scattering in the one-mode approximation is given by Sl]l{q,uj), a sum over inde- 
pendent normal-mode cross sections. Sl]l{q,uj) provides the natural width of the Brillouin 
peak a. For liquid sodium at melt. g,ves a highly aeeurate aceouut of the location 

of the Brillouin peak 51], and as we have seen in the present study, gives an accurate account 
of the Brillouin peak SbTQdb clS well. 

As shown in Sec. II, transits contribute to Fiig{q,t) in two ways. Transit-induced jumps 
in the atomic equilibrium positions and displacements cause decorrelation among the terms 
in Eq. (fTT|l . hence transits enhance the decay of time correlations. Transits also provide an 
additional source of inelastic scattering, hence increase the inelastic cross section. These 
effects are modeled by the strength parameter C{q) and the relaxation function e""^*-'^''* in 
Eq. (fT^ for FR{q,t), and by the relaxation function e""^'-^-'* in Eq. ()14j) for FB{q,t). The 
model so constructed is a generalization of Zwanzig's model for the velocity autocorrelation 

n 

function 35]. The model expressions for Siiq{q,uj) are given in Eqs. ()15II17|1 . 

As shown in Figs. 2-8, the model can be made to fit MD calculations of S{q, u) extremely 
well, almost within computational errors, except for the two largest q values. The small 
inadequacy of the model in the vicinity of the Brillouin peak in Figs. 7 and 8 is apparently 
due to multimode scattering, present in the MD calculation but not in the vibrational theory 
evaluated here. Properties of the fitting parameters are as follows. In the Rayleigh peak 
contribution, C{q) is greater than one (Table I), indicating the presence of inelastic transit 
scattering in the MD calculations. The relaxation rate ai{q) is close to the mean transit 
rate of 2.5 ps^^ as expected (Fig. 9), but the figure suggests that ai(g) approaches zero as 
q ^ 0. The Brillouin peak relaxation rate a2{q) appears to vanish at small q (Fig. 9), which 
would mean that the Brillouin peak width in the liquid is the natural width at small q. 

Finally, while mode coupling theory and V-T theory each provide a physically based model 
capable of accurately fitting S{q, uj) for the liquid at melt, the two models work with entirely 
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different projections of the underlying liquid motion. V-T theory is the more universal, in 
that it applies the same motional constituents to all statistical averages, equilibrium and 
nonequilibrium alike. 
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